function malkov_feature=markov_d(dctresult,T)
    [M,N] = size(dctresult);
    arr = dctresult(1:M-1,1:N-1)-dctresult(2:M,2:N);
    arr(arr>T) = T;
    arr(arr<-T) = -T;
    [row, col] = size(arr);
    m = zeros(2*T+1,2*T+1);
    for i = 1:row-1
        for j = 1:col-1
            v1 = arr(i,j) + T + 1;
            v2 = arr(i+1,j+1) + T + 1;
            m(v2, v1) = m(v2, v1)+1;
        end
    end
    s = sum(m,1);
    for j = 1:2*T+1
        if s(j) > 0
            m(:,j) = m(:,j) / s(j);
        end
    end
    malkov_feature = reshape(m,1,(2*T+1)*(2*T+1));
end
